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ABSTRACT 

We present results of X-ray observations of the Coma cluster with multiple instruments over a broad energy 
band. Using the data from INTEGRAL, RXTE and ROSAT observatories, we find that the Coma spectrum in the 
0.5 - 107 keV energy band can be well approximated by a thermal plasma emission model with a temperature 
of T = 8.2 keV. INTEGRAL was used to image the cluster emission in the hard energy band. The cluster is only 
marginally detectable (~ 1 .6 <x) in the 44-107 keV energy band; however, the raw flux in this band is consistent 
with the previous results from Beppo-SAX and RXTE observatories. We can exclude with high significance that 
the hard-band flux reported by Beppo-SAX and RXTE could be produced by a single point source. The 20 - 80 
keV flux of a possible non-thermal component in the cluster spectrum is (6.0 + 8.8) x 10~ 12 ergs cirT 2 s _1 . It 
is unlikely that the IC scattering of CMB photons is able to produce hard X-ray flux at these levels, unless the 
magnetic field strength is as low as 0.2 pG. The latter value can be considered as a lower limit on the field 
strength in Coma. We also present a temperature map of the central part of the cluster, which shows significant 
variations and in particular, a hot, ~ 11.5 keV, region in the extension towards the subcluster infalling from the 
South- West. 

Subject headings: clusters of galaxies: general — clusters of galaxies: individual(Coma) 



1. INTRODUCTION 

Hot intracluster gas should have an admixture of non- 
thermal components, including relativistic electrons, mag- 
netic fields, relativistic protons, and supra-thermal electrons. 
Some of these components are expected on the theoretical 
grounds while the presence of relativistic electrons and siz- 
able magnetic fields is clear from observations of cluster ra- 
dio halos (see e.g., a recent review by [Ferrari et al. 2008} . 
Often, the very existence of non-thermal components in the 
ICM poses inter esting problems wh ich have been studied in- 
tensively (see Reph aeli et al.|20 08 for a recent review). The 
relativistic electron population in clusters can be studied in 
the hard X-ray band through the high-energy spectral compo- 
nents observed on top of the thermal bremsstrahlung spectra. 
If such components are present at the intensity levels accessi- 
ble with the current instrumentation, their likely origin is in- 
verse Compton (IC) scattering of the cosmic microwave back- 
ground photons; the models in which high-energy X-rays are 
produced by bremsstrahlung of the supra-thermal electrons 
are less attractive (e.g., Petrosian, 2001; Petrosian, Bykov 
& Rephaeli, 2008; Petrosian & Eats, 2008) because subrel- 
ativistic electrons are subject to strong Coulomb losses (but 
see Dogiel et al. ( 2007) >). At some level, an excess around 
20-30 keV with respect to a single-temperature fit can be at- 
tributed to the presence of localized patches of hot gas with 
temperatures substantially exceeding the cluster average. The 
amplitude and the shape of such high energy excess is an im- 
portant diagnostic for the non-isothermality of the gas. 

A classic application of the hard X-ray spectral observa- 
tions is the measurement of the bulk magnetic field strength 
via comparison of the inverse Compton and radio synchrotron 
emission (e.g., Felten & Morrison 1966, Tucker et al., 1973, 
Rephaeli 1979). A prime object for such studies is the Coma 
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cluste r, which posses ses a brig ht, well-s t udied radio halo (see, 
e.g., Thierbach, Klein, & Wielebinski |(|2003|l). Since Coma 



is also the nearest rich cluster, it was observed with virtually 
every X-ray observatory flown. 

Detection of the hard X-ray component in excess over the 
extrapolation of the thermal spectrum in Coma was first re- 
ported from the RXTE and Beppo-SAX observations ( |Rephaeli"] 
|et al.| [l999; Fusco-Femia no et al.||1999|l, later confirmed by 
more extensive analyses (Rephaeli & Gruber 2002} |Fusco-| 
|Femiano et al. 20041. The statistical significance of these 



detections remains not very high (e.g., Fusco-Femiano et al. 
2004, reported a ~ 4.8cr detection), and the reported detec 



tions are subject to criticism (e.g., Ros setti & Molendi|2004| 
report only a ~ 2<x significance of the hard component from 
independent reanalysis of the data). Even if a hard compo- 
nent exists in the Coma spectrum, we cannot exclude, on the 
basis of the Beppo-SAX of RXTE data, a possibility that it is 
produced by, e.g., a strongly absorbed AGN. Obviously, the 
situation can be improved through observations of Coma with 
an imaging hard X-ray telescope. 



Telescopes of the INTEGRAL observatory (Winkler et al. 
|2003| l offer a unique combination of good sensitivity and an- 
gular resolution in hard X-rays. The total INTEGRAL expo- 
sure of the Coma cluster is almost 10 6 sec, leading to the sen- 
sitivity for an extended source at E ~ 50 keV which is com- 
parable to that of the Beppo-SAX and RXTE data. Coma is 
detected by INTEGRAL with high significance at E < 30 keV, 
where the emission is dominated by the thermal component; 
the analysis in t his energy band has been reported by|Renaud| 
|etaLl (T2006 ) and Eck ert et al.| ( |2W) ; [Renaud et al.| ( |2006| l also 
used the first 500 ksec of the Coma data to put limits on the 
hard X-ray component. In this work, we present a systematic 
analysis of the hard X-ray component using the full INTE- 
GRAL exposure in combination with the data from ROSAT 
and RXTE. The combination of the data from these satel- 
lites provides an accurate, self-consistent measurement of the 
broad band thermal spectrum of Coma, which is critical for 
detection of non-thermal components at or below 50 keV. We 
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also present a temperature map of the central region of Coma 
obtained from the ratio of INTEGRAL and ROSAT bright- 
nesses. 

All distance-dependent quantities are reported assuming 
that Coma is at d = 90.5 Mpc (corresponding to z — 0.023 
and H = 70 km s _I Mpc" 1 ). 

2. OBSERVATIONS AND DATA ANALYSIS 

We use X-ray observations of the Coma cluster performed 
by ROSAT, RXTE, and INTEGRAL. These observatories com- 
bined cover a very broad energy range, 0.5-100 keV, and thus 
provide a very high-quality broad band spectrum. Since Coma 
is a bright X-ray source, the spectral measurements with each 
instrument are relatively straightforward. The main compli- 
cation is adjusting the observed spectra to account for the dif- 
ferent spatial response: ROSAT is a direct-imaging telescope; 
INTEGRAL is a coded-mask imager which results in a rather 
broad Point Spread Function (PSF) in reconstructed images; 
RXTE is a collimator with a beam pattern which is slightly 
smaller than the cluster size. Also, the statistical quality of the 
Coma data in the ROSAT and RXTE energy bands is very high 
and so uncertainties are dominated by the instrument cross- 
calibration. These issues are discussed below for each tele- 
scope individually. 

2.1. INTEGRAL 

Coma was targeted by INTEGRAL in several sets of obser- 
vations with a total exposure near 106 seconds (revolutions 
36, 71, 72, 274, 275, 317, 318, 319, 324, 325) We con- 
centrated here on the data from ISGRI detect or (|Lebrun et| 
aLl [2003] i of the IBIS telescope ( |Ubertini et al. 2003 ). This 
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telescope provides a wide field of view (28° x 28°) and an 
angular resolution of 12'. The telescope is sensitive over ap- 
proximately 17 - 200 keV energy band. Some fraction of 
the IBIS data was affected by Solar flares and was discarded 
from the analysis. The clean set of IBIS data was initially 
reduced in individual pointings (so called science windows), 
each with a typical duration ~ 2 ksec. The reconstructed im- 
ages were combined into a single mosaic which was analyzed 
further. The total dead-time corrected combined exposure was 
990 ksec. 

The IBIS/ISGRI image reconstruction algorithm was dis- 
cussed previously in |Krivonos et al. (2005 2007), and we re- 
fer the reader to these works. A crucial prerequisite for this 

— and any other coded mask image reconstruction algorithm 

— is the ability to accurately predict the detector background 
image in the absence of any sources. Our image reconstruc- 
tion software uses standard calibration tables (OSA 6.0 4 ) were 
used to correct the event energies depending of their rise-time. 
We applied, however, additional corrections associated with 



the secular change of gain (see discussion in |Tsygankov et alT 
|2006) . 

Angular size of the Coma cluster is significantly larger 
than the IBIS PSF. This is apparent from the comparison of 
the "growth function" (flux integrated within an aperture) for 
Coma and a point source (Crab) shown in Fig[T] While es- 
sentially all of Crab flux is concentrated within 12' (three 
4' pixels in the reconstructed image), for Coma the curve 
continues to grow up to a 30' radius. In fact, the observed 
growth curve is consistent with that for a /3-model (Cava- 
liere & Fusco-Femiano 1976) fit to the ROSAT image, S{r) = 
S (l+r2/r c 2)-^ +0 \ with/3 = 0.741 and r c = 10.7'. Because 
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Fig. 1 . — The "growth functions" (flux integrated within a given radius in 
the reconstructed image) of the Coma cluster (solid circles) and Crab nebula 
(open circles) in the 17-22 keV energy band. The solid line shows a growth 
curve expected for a yS-model profile with fi = 0.741 and r c = 10.7' (parame- 
ters for the /3-model in the ROSAT band). The Crab flux was rescaled to the 
Coma level for comparison. 

of such a good agreement, the total Coma flux in the INTE- 
GRAL energy bands can be obtained by fitting a normalization 
of the /?-model instead of directly integrating the flux in each 
pixel within a wide aperture. The advantage of this method is 
that it provides a higher statistical accuracy since each pixel 
is added with an optimal weight. Direct integration, however, 
is less model-dependent since the hard-band thermal — and 
especially, the inverse Compton — emission in general does 
not have to follow the distribution of the surface brightness in 
the ROSAT band. Given the pros and cons, our choice is to 
use the /J-model fluxes throughout but always check that they 
are consistent with the direct integration. 

The limiting factor for the INTEGRAL observations of the 
Coma cluster is the level of background fluctuations. In addi- 
tion to purely statistical fluctuations, the background can leave 
systematic variations in the reconstructed image if subtracted 
incompletely from the raw detector image. Every care was 
taken, therefore, to ensure that the background was subtracted 
correctly. Background templates were taken from observa- 
tions of "empty" fields closest in time to each Coma pointing 
or group of pointings. With this approach, we correctly take 
into account all possible temporal variations of ISGRI back- 
ground. As a final check, we checked the statistics of fluctua- 
ti ons in the recons tructed image (similar to the analysis in §3 
of |Krivonos et al.| ( [20 07)), and also fluxes within 40' circles at 
different off-axis locations. In both cases, the variations had 
the mean consistent with zero and the dispersion equal to that 
expected for the purely Poisson fluctuations. We conclude, 
that the accuracy of our background modeling has reached its 
fundamental statistical limit (for the exposure time accumu- 
lated during Coma observations). 

To convert ISGRI counts to physical flux units 
(erg s _1 cm -2 ), we used the calibration observations of the 
Crab nebula which are regularly performing by INTEGRAL. 
Crab observations were reduced with the same software 
setup we used for Coma and the counts-to-flux conversion 
coefficients were determined assuming the "conventional" 
spectral parameters for Crab, / = 9.7 phot s _1 cirT 2 keV~' 
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follows the ROSAT image, the effective flux fraction within 
the PCA FOV can be computed as 
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Fig. 2. — Broadband spectrum of the Coma cluster from INTEGRAL 
(filled circles), RXTE (open squares), and ROSAT (open circles). The best-fit 
MEKAL model is shown by solid line. The contribution of non-thermal com- 
ponents corr esponding to the Rephaeli & Gruber 1 2002 1 and |Fusco-Femiano| 
|et al. |{200 4 1 detections is shown by dotted and dashed lines, respectively. A 
2<r upper limit for a point source emission is shown by a dashed-dotted line. 



( |Toor&Seward|1994| ). 

After experimenting with energy bands for the Coma anal- 
ysis, we chose to extract fluxes in the 17 - 22, 22 - 28.5, 
28.5 - 44, and 44 - 107 keV bands. ISGRI efficiency quickly 
drops below 17 keV setting a natural lower limit for the energy 
band. The width of the two lower channels, factor of 1 .3 in en- 
ergy, was chosen to ensure high statistical significance of the 
Coma detection (> 7<x). The flux in the 44 - 107 keV channel 
should be dominated by non-thermal components. The third 
channel, 28.5 - 44 keV fills the gap. 

Source fluxes were extracted within the 60' radius. This 
aperture, one the one hand, contains almost all of the INTE- 
GRAL X-ray flux (Fig|T]l, and on the other hand matches the 
size of the RXTE field of view, making flux comparisons eas- 
ier. ROSAT pointings fully cover this region. 

2.2. RXTE 

Coma was observed by RXTE for ~ 90 ksec (Obs.ID 10368) 
and the data from its PCA spectrometer can be used to mea- 
sure the cluster spectrum in the 3-20 keV energy band. 
The data reduction was done with standard programs of the 
LHEASOFT package (v6.0). To increase the sensitivity and 
reduce systematic uncertainties, we used only data of the first 
layers of PCA detectors. The background was based on the 
"L7J240CM" model (this includes both the particle-induced 
detector background and the all-sky average of the Cosmic 
X-ray Background). 

One of the complications in the RXTE analysis is to cor- 
rectly compute the flux fraction outside the PCA collimator 
FOV. The PCA beam pattern can be modeled as a convolu- 
tion of a near-conical response of the individual collimator, 
R(x, y), with a Gaussian (<x = 6') which corresponds to mis- 
alignment of individual collimators and spacecraft pointing 
drift during the observation ( |Jahoda et al.| [2006). R(x,y) is 
tabulated in the calibration file pcacol. Assuming that in the 
PCA energy band, the Coma surface brightness distribution 



/ S (x, y) X (R(x, y) <g> G(x, y)) dxdy 
J S (x,y) dxdy 



0.763, 



(1) 



where S(x,y) is the ROSAT image and G(x,y) is a Gaussian. 
To obtain total Coma fluxes, the observed RXTE count rates 
should be divided by this factor. 

In addition to the FOV correction, we need to ensure good 
cross-calibration of the RXTE/PCA with ROSAT and INTE- 
GRAL. To this end, we note that the best-fit parameters of 
the Crab spectrum obtained with the standard PCA response 
matrix differ significantly from the "conventional" values: 
2.09 + 0.04 and 11.6 ± 0.4 phot s _1 cnr 2 keV" 1 instead of 
2.1 and 9.7 phot s _1 crrT 2 keV -1 , respectively, indicating a 



possible problem with the PCA absolute calibration ( Revnivt- 
sev et al. 2003). An ad-hoc correction factor, 



f corr (E) = 0.836£ 



-o.oi 



(2) 



applied to the observed PCA spectra brings the Crab results 
into agreement with the conventional average spectrum ( Toor 



& Seward 1994). We apply this correction factor also to the 
Coma data. We note that this eliminates a 20% difference in 
flux when the PCA spectrum is extrapolated to the ROSAT 
band; the INTEGRAL and PCA fluxes near 20 keV are also 
in a good agreement after this correction is applied (Fig[2]>. 
Since we effectively base the PCA response calibration on the 
absolute Crab spectrum, we need to assign systematic errors 
to account for the uncertainties i n the latter. The unce rtainty 
of the Crab spectral index, +0.04 ( Revnivtsev et al.|2003 1 can 
be approximated if we bin the PCA data into five wide en- 
ergy channels and assign a 7% uncertainty to the flux in each 
channel. 

2.3. ROSAT 

The ROSAT PSPC pointed observations of Coma were re- 
duced as described in Vikhl inin et al.| ( |1999| l. The reduc- 
tion pipeline was based on S. Snowden's software (Snowden 
|et al.|[T994] >. This software eliminates periods of high parti- 
cle and scattered solar backgrounds as well as those intervals 
when the detector may be unstable. Exposure maps in several 
energy bands are then created using detector maps obtained 
during the ROSAT All-Sky Survey. The exposure maps in- 
clude vignetting and all detector artifacts. The unvignetted 
particle background is estimated and subtracted from the data 
even though the PSPC particle background is low compared 
to the cosmic X-ray background. The scattered solar X-ray 
background also should be subtracted separately, because, de- 
pending on the viewing angle, it can introduce a constant 
background gradient across the image. Most of Solar X-rays 
were eliminated by simply excluding time intervals when this 
emission was high, but the remaining contribution was also 
modeled and subtracted. The remaining background was esti- 
mated by extracting the radial surface brightness profile from 
the merged data, and fitting it to the (5 + const model at large 
radii. The background-subtracted images are suitable for di- 
rect extraction of fluxes in the energy bands of interest. Our 
imaging analysis uses the 0.5-2 keV band images. 



3. RESULTS 

3.1. Spectral analysis 

The combined energy spectrum of Coma from INTEGRAL 
(17-107 keV), RXTE (3 - 20 keV), and ROSAT (0.5 - 2 keV) 
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Fig. 3. — Photon spectra of Coma, derived in this work (symb ols with 
the same meaning as in Fig[2j and from Beppo-SAX (open triangles, |Fusco-| 
|Femiano et al. 20041. We dianot adjust the normalization of the Beppo-SAX 
spectrum to account for the cluster flux outside the collimator FOV. Note the 
overall mismatch in spectral slope between our RXTE and INTEGRAL mea- 
surements in the 15-50 keV energy band. This has implications for the de- 
composition of the observed data into thermal and non-thermal components. 

is shown in Figj2] The spectrum can be well fit with the ther- 
mal plasma emission (MEKAL) with the temperature T — 
8.2 ± 0.2 keV, abundance fixed at 0.250 Solar ( |Arnaud etaL] 



2001), and Galactic absorption N H = 9 x 10 ly cm 2 . The 
best-fit approximation is shown in Fig|2]by a solid line. This 
model gives a flux in the INTEGRAL band, 17-60 keV, 
of f x = 5.6 x 10" 11 erg cm" 2 s" 1 and luminosity L x - 
7.1 x 10 43 ergs s _1 . The corresponding flux and luminosity in 
the 2-10 keV energy band are f x = 4.2 x 10" 10 erg cm" 2 s" 1 
and L x = 5.1 x 10 44 ergs s _1 . The flux and luminosity un- 
certainties in these bands are dominated by systematic errors 
which we estimate to be of order 5%. 

INTEGRAL fluxes at E < 44 keV are in full agreement with 
the thermal spectrum. The thermal component (including pos- 
sible contribution from hotter regions within the cluster, see 
below) should make a negligible contribution at higher ener- 
gies. The INTEGRAL-measmed flux in the 44-107 keV band 
is (1.8±l.l)xl0 -11 ergs cm" 2 s . If real, this emission cannot 
be attributed to the high energy tail of the thermal spectrum 
and should instead represent a non-thermal component. We 
note, however, that the statistical significance is low (1.6 cr) 
and INTEGRAL statistical uncertainties are comparable to, or 
higher than, those in the earlier RXTE or Beppo-SAX measure- 
ments (|Rephaeli et al.|19"99"l|Rephaeli & Gruber|2002| ~ 



Femiano et al. 
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1999 2004J7 We also tried to repeat the 



Fusco- 

( 2004) procedure by fixing the photon index of 



Fusco- 



a non-thermal component at T = 2 and fitting the total spec- 
trum with a thermal plus power law model. This gives a power 
law flux of (6.0 ± 8.9) x 1 0~ 12 ergs cm" 2 s" 1 in the 20- 80 keV 
energy band (same as in Fusco-Femiano et al.|2004|l. This is 
a factor of 2.5 lower than the flux reported by|Fusco-Femiano 



et al. ( 2004 1 in the same energy band even though we INTE- 
GRAL and Beppo-SAX appear to measure the same flux above 
~ 50 keV (Fig(3]). The source of discrepancy seems to be 
instead in the fluxes measured at the lower boundary of the 



Beppo-SAX band. The Beppo-SAX data points below 20 keV 
seem to be a factor of 1 .5-2 below our RXTE and INTEGRAL 
measurements, and hence they imply a lower thermal compo- 
nent flux. These discrepancies underscore the importance of 
using a broad-band spectra for detection of non-thermal com- 
ponents from hot clusters at around E = 50 keV. 

Where INTEGRAL can do significantly better than the pre- 
vious observatories is to check whether this emission can be 
attributed to single point source within the cluster. We do 
not detect any significant point sources in the hard-band im- 
age (see below) and a 2cr upper limit on the point source flux 
is 7.6 x 10" 12 ergs cm" 2 s" 1 in the 44-107 keV band. This 
is below the fluxes of non-thermal components reported by 
Fusco-Femiano et al.| ( 2004} and Rephaeli & Gruber ( 2002 1 



and so we can exclude a possibility that these detections can 
be attributed to a single persistent AGN. 

3.2. Imaging analysis 

To study the spatial structure of the cluster emission, we 
used images in the 17 - 28.5 keV ("soft") and 44 - 107 keV 
("hard") INTEGRAL bands. The first band is a combination of 
two spectral channels where the cluster emission was detected 
with high significance; the second band is where the putative 
non-thermal component should dominate the thermal plasma 
emission. The images are shown in Fig |4] In the soft band, 
Coma is clearly an extended source (see Eckert et al. 20071 for 



detailed modeling of the INTEGRAL image); a point source 
would be confined to a ~ 3 x 3 pixels square in these images. 

The contours in Fig|4]show the ROSAT surface brightness 
levels in the 0.5-2 keV band. The INTEGRAL soft-band im- 
age shows an elongation towards the infalling subcluster and 
a small offset relative to the ROSAT surface brightness peak 
(see also |Eckert et al.|2007) . The off-set is ~ 4.3' to the West 
(FigQ; it is small but significant (e.g., the locations of the two 
AGNs detected in the Coma field, NGC 4151 and NGC 4388 
are within 0.2' of their optical positions). As we discuss be- 
low, the offset between INTEGRAL and ROSAT images most 
likely reflects the temperature variations within the cluster. 

We now discuss the hard-band INTEGRAL image (44- 
107 keV). As was discussed above, the total flux within 60' 
was detected in this band with a 1 .6 cr significance. Can this 
emission be associated with a small number of point sources 
or with an extended component centered on the cluster? The 
raw image (Pig J4j?) does not show any significant structures. 
The brightest spots correspond to 1 .5 - 2 cr significance. In 
particular, we can exclude a possibility that the flux measured 
by INTEGRAL in this band can be attributed to a single point 
source; a 2 cr upper limit is ~ 7.6 x 10" 12 ergs cm" 2 s . Fig- 
ure [5] shows the hard-band image smoothed with the ROSAT 
y6-model to illustrate the distribution of a more extended com- 
ponent. Again, there are no statistically significant features in 
the smoothed image. A local maximum near the Coma center 
is offset from the cluster location and in fact is not the bright- 
est feature in the field. Therefore, we cannot conclude that the 
hard emission is spatially coincident with the Coma cluster; it 
is more consistent with being a statistical fluctuation. 

3.3. Temperature Map 

INTEGRAL'S ability to reconstruct X-ray images in the 20- 
30 keV band opens a unique opportunity to compute the 
temperature map of Coma simply from the ratio of INTE- 
GRAL and ROSAT surface brightnesses. The INTEGRAL 
energy band is above the exponential cutoff in the thermal 
T ~ 10 keV spectrum, thus the resulting temperature map 
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Fig. 5. — Same as in FigHri, but converted with the beta-model. The gray 
scale is different than in Figpfk 



Fig. 4. — INTEGRAL images (~ 3° X 3°) of the Coma region obtained in 
17 - 28.5 keV (a) and 44 - 107 keV (b) energy ranges. Contours represent 
the ROSAT intensity levels in the 0.5 - 2 keV energy band. Circles indicate 
an r = 60' regions around the cluster. 



is very robust and relatively insensitive to the calibration un- 
certainties in either instrument. 

The drawback is a relatively poor PSF of INTEGRAL/IBIS 
which results in a poor spatial resolution in the resulting tem- 
perature map. The PSF angular size is comparable to the clus- 
ter core-radius and so the ROSAT -to-INTEGRAL ratio should 
take into account the redistribution of flux from the center 
to large radii. To improve statistics, we smoothed the IN- 
TEGRAL image (Fig|4^) with a Gaussian with cr = 5' that 
approximates the IBIS PSF in the mosaic images contain- 
ing observations with different rotation angles (Krivonos et al. 
2007) . The smoothed image has an effective PSF that can be 
approximated by a cr = 7.1' Gaussian. To match the resolu- 
tion of soft- and hard-band data, we smoothed the ROSAT im- 
age with the same Gaussian before computing the ratio map. 




Fig. 6. — The temperature map of the central part of the Coma cluster. Tem- 
peratures of 7.5-8.5 keV, 9.5-10.5 keV and more than 1 1 keV correspond to 
brown, yellow and white colors, respectively. Contours represent the ROSAT 
intensity levels. 

The obtained INTEGRAL-Xo-ROSAT ratio map was con- 
verted to the temperature map using a lookup table of 17-28.5 
to 0.5-2 keV flux ratios for the MEKAL model as a function 
of temperature. As can be expected, this ratio is very sensitive 
to the temperature (varies by a factor of ~ 10 when T varies 
from 5 to 12 keV). An approximately 10 cr detection of Coma 
by INTEGRAL translates into ~ 0.3 keV temperature uncer- 
tainties in the central region. The resulting temperature map 
is shown in Fig|6] The map is restricted to the region where 
the hard X-ray flux is detected with at least 2.5 cr significance. 
Qualitatively, our temperature map is similar to that obtained 
by lArnaud et al.l (1200 lb from the XMM-Newton data. The 
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most notable features are the cold region to the South-East of 
the center coinciding with a filamentary structure in the X-ray 
brightness (Vikhlinin et al. 1997\ and the hot region in the 
direction of the infalling subcluster 5 . 

4. DISCUSSION 

The possibility of using radio (synchrotron) and X-ray (In- 
verse Compton) observations to constrain the strength of the 
magnetic fields in plasma has been extensively discussed in 
application to various astrophysical sources (e.g. Felten & 
Morrison 1966, Tucker et al., 1973, Rephaeli 1979). We 
make similar estimates using the parameters relevant for IN- 
TEGRAL observations of the Coma cluster. For an electron 
with the Lorentz factor y » 1 moving in a uniform magnetic 
field B at a pitch angle 8 the synchrotron emission spectrum 
peaks at the frequency 



0.29- 



3 eBsinO 
An m e c 



-r 



(3) 



(e.g., Ginzburg & Syrovatskii 



electron mass and charge, and t 



1965), where m e , e, c are the 



le speed of light, respectively. 



The synchrotron emission of the Coma halo was observed at 
frequencies ranging from 30 MHz up to 1.4 GHz, and was 
shown to have spectral index a ~ 1.34 (see e.g. Kim et al., 
1990, Deiss et al., 1997 and references therein). At higher 
frequencies the evidence for spectrum steepening has been 
reported (Schlickeiser, Sievers & Thiemann 1987, Thierbach, 
Klein, & Wielebinski 2003). Below we will use the halo flux 



F r = .64 ± 0.035 Jy at v r = 1.4 GHz reported by Dei ss et al. 
( 19971. From eqj3]it follows that emission at this frequency 
is provided by electrons with y r ~ 10 4 if the field is order of 
10 pG. We further assume that the surface brightness distri- 
bution at 1.4 GHz can also be approximated by a /3-model; 



S r (x)K[l+(x/r c , r ) 2 ] 



-3B r +0.5 



(4) 



where r c r and B r are the co-radius and beta-parameter of the 
radio surface brightness. According to Deiss et al. ( 1997 1 the 
scale size of the radio halo at 1 .4 GHz is similar to that of the 



X-ray halo parameters derived by |Briel, Henry & Boehringer 
[ ( 1992) > (i.e. r Ctr ~ 10/5, but the radio declines with ra- 
dius more rapidly (i.e. B r > B x = 0.75). |Colafrancesco 



Marchegiani & Perola ( 2005 } fitted the surface brightness 
distribution from Deiss et al.l (11997b with r cr — 23'. 8 and 



B r = 1.47. The volume emissivhy T(r) of the cluster (assum- 
ing spherical symmetry) at 1 .4 GHz is then 



T r (r) oc [l + {rlr c , r f] 



-3B, 



(5) 



We now can estimate the hard X-ray flux due to Inverse 
Compton scattering of CMB photons by relativistic electrons. 
The electron with the Lorentz factor y » 1 will up-scatter 
CMB photons to a characteristic frequency 



4 2 4 2 3kT CMB 

Vx ~ -ry V C MB ~ ; , 

3 3 h 



(6) 



where Tqmb ~ 2.7K is the CMB temperature and k is the 
Boltzmann constant. Thus for the 75 keV photon the required 
electron Lorentz factor is yx ~ 9 x 10 3 , i.e. close to the value 



This hot region was also identified by |Eckert et al|j2007} with a similar 
technique (ratio of INTEGRAL and XMM-Newton fluxes); their temperature 
for this region is however higher than ours (~ 12 keV vs ~ 9.7 keV at the dis- 
tance of ~ 14' from the center to south-west) probably because of neglecting 
the INTEGRAL PSF effects 
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Fig. 7. — Bottom panel: Lorentz factor of electrons emitting synchrotron 
radiation at 1.4 GHz (y r , solid line) as a function of magnetic field B and 
the Lorentz factor of electrons producing 75 keV photons by up-scattering 
CMB radiation (yx, dotted line). Middle panel: The life time of electrons 
having Lorentz factors y r yx (solid and dotted lines respectively). The top 
panel shows the comparison of the spectral intensity at 75 keV (vF v , in units 
of ergs cm -2 s , horizontal dotted line) corresponding to the observed INTE- 
GRAL 44-107 flux, (1 .8 ± 1 . 1) X KT 1 1 ergs cnT 2 s" 1 , with the expected level 
of IC emission. The model IC fluxes were computed using eq[7jfor the halo 
radio flux F r = 0.64 Jy at 1 .4 GHz and uniform magnetic field~B. The short 
and long dashed lines show the same IC level, but scaled by factors of 2.8 
and 1 1 to illustrate the changes introduced by the assumption that magnetic 
fields energy density declines with radius according to eq[8]with the parame- 
ter a = 0.5 and a = 1 respectively. In this case the value ofo shown along X 
axis is the value of magnetic field Bo at the center of the cluster. 

needed to produce radio emission at 1.4 GHz for plausible 
values of magnetic fields (see Figj7] bottom panel). The ex- 
pected IC flux due to CMB photons up-scattering prod uced by 
the electrons responsible for synchrotron emission is (Felten 
|& Morrison [T9661 ): 

F x = F r xC(p)(^j T^ B B- a -\ (7) 

where C{p) is a function of the power law slope of the elec- 
tron energy spectrum, p, which is related to the synchrotron 
spectral index as p — 2a + 1 . 

Assuming that B is constant across the cluster, one can cal- 
culate expected IC flux for a given B and radio flux F r using 
eq|7] as is done in the top panel of Fig 7] The intersection 
of the predicted (solid line) and observed (dotted line) curves 
gives the value of the magnetic field ~ 0.1-0.2 pG needed to 
explain radio and hard X-ray fluxes as the synchrotron and IC 
emission respectively produced by the same electron popula- 
tion. 

The assumption of constant B across the cluster probably 
is too simplistic, especially given that energy density of the 
ICM varies strongly from the center to outskirts (see, e.g., 
Dolag, Bartelmann & Lesch, 2002). The spatial variations of 
the magnetic field would affect the relation of the synchrotron 
and IC fluxes (Rephaeli 1979; Goldshmidt & Rephaeli 1993; 
Brunetti et al. 2001; Newman et al. 2002; Colafrancesco et 
al. 2005). We estimate the effects of magnetic field variations 
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Fig. 8. — Variations of IC flux on the different behavior of the magnetic field 
with radius. The parameterization is done in the form S 2 (r) = ■ See 

text for details. 

using a simple parameterization, 

B\r) = Bl(p{r)lp ) a , (8) 

where p(r) is the thermal gas density, which scales approxi- 
mately as the ICM pressure if the gas density variations are 
much stronger than those of temperature over the region of 
interest. Adopting the ROSAT /?-model for the gas density 
distribution, we can write: 



B(r) = B [l + (x/r c , x ) 2 ] 



-J aft 



(9) 



where Bo is the magnetic field at the center of the cluster. In 
this parameterization a = corresponds to the case of a con- 
stant magnetic field, while a - I implies that energy density 
of the magnetic field scales as the ICM pressure. Compared 
to the constant B case (i.e. eq 7]) the ratio of the IC and syn- 
chrotron fluxes (within given distance from the center) has to 
be multiplied by an additional factor f{a): 



f(a) = 



/ r r (r)r 2 dr 



(10) 



where T r (r) is given by eq|5] where the integration is in the 
range r = 0- 60'. The computed factor f(a) is shown in FigjH] 
(in these calculations we assumed the X-ray surface bright- 
ness parameters r cr = r cx = 10'.68, /3,. = B x — 0.741). 

The total IC flux (up to infinity) may actually diverge for 
large a since the drop of the B at large distances from the cen- 
ter has to be "compensated" by energy density of relativistic 
particles to maintain the assumed TE(r). This causes strong 
increase of IC emission. Formally for convergence (up to in- 
finity) one needs: 



a < 



2(2B r - 1) 
B x (a+l)' 



(11) 



which translates to a < 0.55 for our parameters. However, 
since we integrate the hard X-ray flux only within a 1° radius, 



we can use wider range of a. It follows from FigjH] that the 
change in a from to 1 changes the ratio f{a) by a factor of 
~11. This translates into the uncertainty in the estimated Bo 
by a factor of f(a) l ^ l+a ^ ~ 2.8. An additional uncertainty is 
introduced by our determination of the hard X-ray flux which 
was done by convolution of the INTEGRAL image with the 
surface brightness model observed at E = 1 keV by ROSAT. 
If, for example, the true surface brightness in the hard band 
has the same core-radius but /3 is varied in the range 0.3-1.5, 
this would change our derived values of Bo by +30%. 

Although the modeling uncertainties described above are 
substantial, they do not change the conclusion that it is un- 
likely that the IC emission in the 44-107 keV band can reach 
the flux levels comparable to the ISGRI//JVr£'G/ML sensitiv- 
ity, unless the magnetic field is as low as a fewxO.l /iG. For 
these low values of the magnetic field the life time is shorter 
for electrons responsible for the radio synchrotron emission 
(Pigj7j and thus the IC spectrum of the electron population 
has to be close to a power law, as long as a power law spec- 
trum is observed in the radio. 

If we treat the observed hard X-ray flux as a detection, the 
derived B ~ 0.1 - 0.2 juG is an order of magnitude below the 
values derived from Faraday rotation B ~ 1.7 ji/G (e.g., Kim 
et al.| 1990 >. Various ways to bring these estimates into agree- 
ment have been discussed (see e.g. Goldshmidt & Rephaeli 
1993, Petrosian 2001, Brunetti et al., 2001, Newman et al. 
2002). The most obvious factor affecting these estimates is 
the spatial variations of the field strength. Note that the Fara- 
day rotation signal depends on the line of sight integral of the 
product of the electron density and the parallel component of 
the field. The value of the field derived from the compari- 
son of the cluster IC and radio fluxes depends instead on the 
volume averaged non-linear function of the field (see eq[T0|). 
Thus different assumptions on the non-uniformity of the mag- 
netic field will affect differently the estimates of B with these 
two methods and the above mentioned discrepancy can be at 
least partly removed. 

However, given the low significance of the signal, we 
should treat the INTEGRAL results as an upper limit on the 
IC flux. Our 2<x upper limit in the 44 - 107 keV band is 
perfectly consistent with the magnetic fields stronger than 
fewxO.l fiG. For even larger fields, the expected IC signal 
drops as * B~ a =s B 2 M and will be below the reach of even 
planned future experiments. For example, for B - 10 fiG, the 
expected IC flux at E — 75 keV is a factor of ~ 10 4 below the 
INTEGRAL sensitivity, even without account for any possible 
electron aging (see Fig|7]l. 

5. SUMMARY 

The summary of main results from the deep observations 
of the Coma cluster with INTEGRAL and combining these 
observations with results of RXTE and ROSAT observatories 
lead us to the following conclusions: 

1) The total cluster spectrum in the 0.5-50 keV energy band 
is well described by the thermal plasma emission with the 
mean temperature T 8.2 keV; 

2) There are significant temperature variations within the 
cluster and the mean 8.2 keV temperature is the result of 
mixing emission components with T ranging from 7.5 to 
10.5 keV. The temperature near the ROSAT surface brightness 
peak is ~ 8.5 keV. 

3) We do not detect a significant excess over the thermal 
spectrum at high energies. The upper limit on the non- thermal 
flux is, however, consistent with the previous detections re- 
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ported on the basis of Beppo-SAX and RXTE observations. 

4) It is unlikely that IC emission in hard X-rays can reach 
the flux level comparable with the INTEGRAL/lSGRl sensi- 
tivity for the magnetic fields strongerr than fewxO.l pG. A 2cr 
upper limit of the hard X-ray flux from INTEGRAL is consis- 
tent with B < 1 pG. 

This work is based on observations with INTEGRAL, and 
made use of the INTEGRAL Science Data Center (Versoix, 
Switzerland), Russian Science Data Center of INTEGRAL 



(Moscow, Russia) and the High Energy Astrophysics Sci- 
ence Archive Research Center Online Service, provided by 
the NASA/Goddard Space Flight Center. The work was sup- 
ported by the Chandra grant G05-6121 A, RFBR grant 07-02- 
01051, DFG grant CH389/3-2, the "Extended objects in the 
Universe" program of the Russian Academy of Sciences, and 
a program of support for leading scientific schools (Project 
NSH-1 100.2006.2). AL thanks CfA for hospitality during the 
course of this research. 
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